UK/12-12 
PUPT-2431 



Quantum Quench Across a Zero Temperature 
Holographic Superfluid Transition 



> 



O Pallab Basu (a Q Diptarka Das (a fl Sumit R. Das (a @ and Tatsuma Nishioka 

> 
O 

(a) Department of Physics and Astronomy, 
University of Kentucky, Lexington, KY 40506, USA 
^Department of Physics, 
Princeton University, Princeton, NJ 08564, USA 

-t— > 

D Abstract 

l — 1 We study quantum quench in a holographic model of a zero temperature insulator- 

superfluid transition. The model is a modification of that of ar Xiv:0911.0962| and 
involves a self-coupled complex scalar field, Einstein gravity with a negative cosmo- 
C ■ logical constant, and Maxwell field with one of the spatial directions compact. In a 

j- — suitable regime of parameters, the scalar field can be treated as a probe field whose 

— backreaction to both the metric and the gauge field can be ignored. We show that 

^ when the chemical potential of the dual field theory lies between two critical values, 

i—i the equilibrium background geometry is a AdS soliton with a constant gauge field, 

while the complex scalar condenses leading to broken symmetry. We then turn on a 
^ time dependent source for the order parameter which interpolates between constant 

values and crosses the order-disorder critical point. In the critical region adiabatic- 
ity breaks down, but for a small rate of change of the source v there is a new small-v 
expansion in fractional powers of v. The resulting critical dynamics is dominated by 
a zero mode of the bulk field. To lowest order in this small-?; expansion, the order 
parameter satisfies a time dependent Landau-Ginsburg equation which has z = 2, 
but non-dissipative. These predictions are verified by explicit numerical solutions 
of the bulk equations of motion. 
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1 Introduction and summary 



Recently there has been several efforts to understand the problem of quantum or thermal 
quench [T[j3] in strongly coupled field theories using the AdS/CFT correspondence |3}j7]. 
This approach has been used to explore two interesting issues. The first relates to the 
question of thermalization. In this problem one typically considers a coupling in the 
hamiltonian which varies appreciably with time over some finite time interval. Starting 
with a nice initial state (e.g. the vacuum) the question is whether the system evolves into 
some steady state and whether this steady state resembles a thermal state in a suitably 
defined sense. In the bulk description a time dependent coupling of the boundary field 
theory is a time dependent boundary condition. For example, with an initial AdS this 
leads to black hole formation under suitable conditions. This is a holographic description 
of thermalization, which has been widely studied over the past several years [8 10 with 
other initial conditions as well. 

Many interesting applications of AdS / CFT duality involve a subset of bulk fields whose 
backreaction to gravity can be ignored, so that they can be treated in a probe approxi- 
mation. One set of examples concern probe branes in AdS which lead to hypermultiplet 
fields in the original dual field theory. Even though the background does not change in the 
leading order, it turns out that thermalization of the hypermultiplet sector is still visible 



this manifests itself in the formation of apparent horizons on the worldvolume 11, 12 



The second issue relates to quench across critical points |1}|3]. Consider for example 
starting in a gapped phase, with a parameter in the Hamiltonian varying slowly compared 
to the initial gap, bringing the system close to a value of the parameter where there would 
be an equilibrium critical point. As one comes close to this critical point, adiabaticity 
is inevitably broken. Kibble and Zurek jT, 13 14 argued that in the critical region the 
dynamics reflects universal features leading to scaling of various quantities. These ar- 
guments are based on rather drastic approximations, and for strongly coupled systems 
there is no theoretical framework analogous to renormalization group which leads to such 
scaling. For two-dimensional theories which are suddenly quenched to a critical point, 
powerful techniques of boundary conformal field theory have been used in |3] to show 
that ratios of relaxation times of one point functions, as well as the length/time scales 
associated with the behavior of two point functions of different operators, are given in 
terms of ratios of their conformal dimensions at the critical point, and hence universal. 

In 15 quench dynamics in the critical region of a finite chemical potential holographic 
critical point was studied in a probe approximation. The "phenomenological" model used 
was that of [i~6] which involves a neutral scalar field with quartic self-coupling with a 
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mass-squared lying in the range —9/4 < m 2 < —3/2 in the background of a charged AdS^ 
black brane. The self coupling is large so that the backreaction of the scalar dynamics 
on the background geometry can be ignored. The background Maxwell field gives rise 
to a nonzero chemical potential in the boundary field theory. In 16 it was shown that 
for low enough temperatures, this system undergoes a critical phase transition at a mass 
m 2 . For m 2 < m 2 the scalar field condenses, in a manner similar to holographic super- 
fluids 17-22 . The critical point at m? = m 2 c is a standard mean field transition at any 
non-zero temperature, and becomes a Berezinski-Kosterlitz-Thouless transition at zero 
temperature, as in several other examples of quantum critical transitions. In [15] the crit- 
ical point was probed by turning on a time dependent source for the dual operator, with 
the mass kept exactly at the critical value, i.e. a time dependent boundary value of one of 
the modes of the bulk scalar. The source asymptotes to constant values at early and late 
times, and crosses the critical point at zero source at some intermediate time. The rate 
of time variation v is slow compared to the initial gap. As expected, adiabaticity fails as 
the equilibrium critical point at vanishing source is approached. However, it was shown 
that for any non-zero temperature and small enough v, the bulk solution in the critical 
region can be expanded in fractional powers of v. To lowest order in this expansion, the 
dynamics is dominated by a single mode - the zero mode of the linearized bulk equation, 
which appears exactly at m 2 = m 2 . The resulting dynamics of this zero mode is in fact a 
dissipative Landau-Ginsburg dynamics with a dynamical critical exponent z — 2, and the 
order parameter was shown to obey Kibble-Zurek type scaling. 



The work of 15 is at finite temperature - the dissipation in this model is of course 
due to the presence of a black hole horizon and is expected at any finite temperature. It is 



interesting to ask what happens at zero temperatures. It turns out that the model of 16 



used in 15 becomes subtle at zero temperature. In this case, there is no conventional 
adiabatic expansion even away from the critical point (though there is a different low 
energy expansion, as in [23]). Furthermore, the susceptibility is finite at the transition, 
indicating there is no zero mode. While it should be possible to examine quantum quench 
in this model by numerical methods, we have not been able to get much analytic insight. 

In this paper we study a different model of a quantum critical point, which is a variation 
of the model of insulator-superconductor transition of 24 . The model of 24 involves a 
charged scalar field minimally coupled to gravity with a negative cosmological constant 
and a Maxwell field. One of the spatial directions is compact with some radius R, and 
in addition one can have a non-zero temperature T and a non-zero chemical potential \i 
corresponding to the boundary value of the Maxwell field. In the absence of the scalar 
field this model has a line of Hawking-Page type first order phase transitions in the T-fi 
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plane which separates an (hot) AdS soliton and a (charged) black brane. Exactly on the 
T = line, the two phases correspond to the AdS soliton with a constant Maxwell scalar 
potential, and an extremal black hole. In 24 it was shown that in the presence of a 
minimally coupled charged scalar, the phase diagram changes. When the charge is large 
the scalar and the gauge fields can be regarded as probe fields which do not affect the 
geometry. Now there is a phase with a trivial scalar and a phase with a scalar condensate. 
In the boundary theory the latter is a superfluid phase. This phase transition persists at 
zero temperature, where it separates an unbroken phase at low chemical potential and a 
broken phase - in both cases the background geometry is the AdS soliton, while the gauge 



field is non-trivial in the superfluid phase. The phase diagram is given in Figure 9 of 24 



The idea now is to probe the dynamics of this insulator-superfluid transition at zero 
temperature by turning on a time dependent source for the operator dual to the charged 
field. So long as the scalar is minimally coupled and the charge q is large, this would 
involve analyzing a coupled set of equations of the scalar field and the gauge field. 

However, it turns out that a slight modification of the model allows us to ignore the 
backreaction of the scalar to the gauge field as well. This involves the introduction of 
a quartic self coupling of the scalar A. Then in the regime A ^> q 2 and A ^> k 2 (where 
K is the gravitational coupling), we can consider the dynamics of the charged scalar in 
isolation. 

In this work we first show that in this regime of the parameters the insulator-superfluid 
transition persists. Concretely, for a sufficiently small negative m 2 , there is a critical value 
of the background chemical potential beyond which a nontrivial static solution for the 
scalar becomes thermodynamically favored. Note that unlike other models of holographic 
superconductors the trivial solution does not become dynamically unstable. Rather the 
non-trivial solution has lower energy. The transition is a standard mean field critical 
transition. The background geometry remains an AdS soliton and the background gauge 
potential remains a constant, which is the chemical potential \x. At the transition, the 
linearized equation has a zero mode solution which is regular both at the boundary and 
at the tip. 

We then turn on a time dependent boundary condition and find that the breakdown 
of adiabaticity for a small rate v is characterized by exponents which are appropriate 
for a dynamical critical exponent z — 2. In a way quite similar to [15| we find that in 
the critical region there is a new small v expansion in fractional powers of v, and the 
dynamics is once again dominated by a zero mode. The real and imaginary parts of the 
zero mode now satisfy a coupled set of Landau-Ginsburg type equation with first order 
time derivatives. However the resulting system is oscillatory rather than dissipative - this 
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is expected since the background geometry has no horizon so that we have is a closed 
system. The order parameter is shown to obey a Kibble-Zurek type scaling. Finally we 
solve the bulk equations numerically and verify the scaling property obtained from the 
above small-f expansion. 

Thermal quench in holographic superfluids with backreaction has been recently studied 
in 25 . This work addresses a different issue - here the quench is applied to the system 
in the ordered phase away from the critical point and the resulting late time relaxation 
of the order parameter is studied. Our emphasis is on probing a possible Kibble-Zurek 
scaling when the quench crosses the critical point. 

In Section 2 we define the model and discuss its equilibrium phases. In Section 3 we 
study quantum quench in this model by turning on a time dependent source, discuss the 
breakdown of adiabaticity and show that the critical region dynamics is dominated by the 
zero mode, leading to scaling behavior. In Section 4 we present the results of a numerical 
solution of the equations, verifying the scaling behavior. In an appendix we discuss a 
Landau-Ginsburg model similar to the critical dynamics of our holographic model. 



2 The model and equilibrium phases 

The "phenomenological" holographic model we consider is a slight variation of the model 



of 24 . The bulk action in (d + 2)-dimensions is 



S = J d d+2 x^j 



i f R+ d(d+iy\ _ i _ i / _ . 2 _ m2|$|2 _ i 4 



(2.1) 

where $ is a complex scalar field and is an abelian gauge field, and the other notations 
are standard. Henceforth we will use L — 1 units. 

One of the spatial directions, which we will denote by 9 will be considered to be 
compact. We will consider the regime 

A > q 2 , A > k 2 . (2.2) 

In this regime the scalar field is a probe field, and its backreaction to both the metric and 
the gauge field can be ignored. 

2.1 The background 

The background metric and the gauge field can be then obtained by solving the Einstein- 
Maxwell equations with the appropriate periodicity condition on 6. It is well known that 



5 



there are two possible solutions. The first is the AdSd+2 soliton, 

d-l 



ds 2 
fsi(r) 



dr^ 



r 2 f s i{r) 

r \d+i 



+ r 2 -dt 2 + dx 2 + r 2 f sl (r)d6 2 



i=l 



with constant parameters \i and r . The periodicity of 9 in this solution is 

6 ~ 6 + 



(d + l)r 



(2.3) 



(2.4) 



while the temperature can be arbitrary. The second solution is a AdSd+2 charged black 
hole 



ds 2 



-r 2 f bh {r)dt 2 



dr 2 



'd-l 



r 2 fbh(r) 



+ r 2 [Y^dx 2 + de 2 



1 



d — 1 / fi 



2d \r 



d-l 



hh{r) = 1- 
A t = fx 

The temperature of this black brane is 



r. y'- 1 d-l / // \ 2 //•. n 
r J 2d \r+ J V r 



(2.5) 



47T 



d+1 - 



(d-l) 2 



2d Vr 



(2.6) 



while the period of 9 is arbitrary. As shown in 24 , this system undergoes a phase 
transition between these two solutions when 



„d+l _ d+l 



1 + 



d — 1 f fi 



2d \r + 



(2.7) 



The AdS soliton is stable when the temperature and the chemical potential are small. At 
T = the transition happens at a critical chemical potential fi C 2 given by 



r {d+ l)(2d) W) 

(d-l)3TX(d+ 1)V2 



(2.J 
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2.2 Scalar condensate 



Consider now the scalar wave equation in the AdS soliton background (2.3). We first 
rescale 

r ->■ — , t ->■ tr , /i ->• — . (2.9) 

In the rest of the paper we will use these rescaled coordinates (i.e., r = 1) and chemical 
potentials. 

For fields which depend only on t and r, the equation of motion is given by 



-Ud t -i^ + ^d r (r^f sl (r)d r ) 



$ - m 2 $ - $|$| 2 = . (2.10) 



In this paper we will consider — ^ d+ ^ < m 2 < _ The asymptotic behavior of the 

solution at the AdS boundary r — > oo is of the standard form 

$(r, t) = J(i) r~ A ~ [1 + 0(l/r 2 )] + A(t)r~ A+ [1 + 0(l/r 2 )] + • • • , (2.11) 

where 

± = <±± ± f^V. ,2,2) 

In "standard quantization" J(t) is the source, while the expectation value of the dual 
operator is given by 

(0} = A(t). (2.13) 

In "alternative quantization" the role of J(t) and A(t) are interchanged. In this mass 
range both A± are positive and both the solutions of the linear equation vanish at the 
boundary. Thus the nonlinear terms in the equation (2.10) are subdominant - which is 
why the leading solution near the boundary is the same as those of the linear equation, 
as written above. 



We need to find time independent solutions of the equation (2.11). Because of gauge 
invariance, we need to specify a gauge to qualify what we mean by time independence. 
For the equilibrium solution we require the solution to be real - this fixes the gauge. Note 
that the tip of the soliton is locally two-dimensional flat space. Therefore we need to 
require the solution to be regular at the tip r = 1. This leads to the following boundary 
condition at r = 1 

$(r)=$ fc + $' ft (r-l) + ... , (2.14) 

where regularity requires 
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To examine the phase structure we need to find time independent solutions with a van- 
ishing source. 

Clearly $ = is always a solution. We have solved the equations numerically and 
found that there is a critical value of the chemical potential [i c \ beyond which there is 
another solution with a non-trivial r dependence which is thermodynamically preferred. 
This means that for // > /i c i, the operator dual to the bulk scalar has a vacuum expectation 
value, i.e., the global U(l) symmetry of the boundary theory is spontaneously broken. 
Although this could happen both in the standard and alternative quantizations, we need 
to check the critical value is less that that of the phase transition between the AdS soliton 
and AdS black hole: /i c i < /i C 2. Otherwise, the scalar condensate phase is not available 
on the AdS soliton. 

Figure [I] shows the behavior of the expectation value (O) for m 2 = —15/4 for standard 
and quantization. We are plotting the condensation with respect to [iq and the phase 
transition happens at fi c \q ~ 1.89, which means the critical chemical potential is very 



small of order 0(l/q) in the probe limit. It follows from (2.8) that /i c i is always much 
smaller than /i c2 ~ 1.86 and there exists a scalar condensate phase on the AdS soliton. 
Similarly for any given m 2 , fi = fi c i is a critical point by letting q be large enough. 



<0> 




Figure 1: The condensations of the scalar operators. 



This transition was first found in 24 for a minimally coupled complex scalar - in this 



case the backreaction to the gauge field cannot be ignored, and the result follows from an 
analysis of the coupled set of equations for the scalar and the gauge field. In this case the 
gauge field introduces non-linearity in the problem which is necessary for condensation of 
the scalar. What we found is that a self-coupling does the same same job. 



8 



2.3 The zero mode at the critical point 



To get some insight into this transition it is useful to write the equation (2.10) as a 
Schrodinger problem. First define a new coordinate 

P(r)=r^TT, (2-16) 

which is the "tortoise coordinate" for the AdS soliton. p(r) is a monotonic function of r 
with the behavior 

p ~ 1/r , r — )■ oo , 

oJr — 1 

p -> p, + ^^, r ^l. (2.17) 
V« + 1 

For example, for d = 3 (asymptotically AdSs spacetime) soliton p* = 1.311. Let us now 
redefine the field by 

_1 fdp^ 1/2 
[r(p 

Then ^(p, t) satisfies the equation 



r 2-d 



\-dl + 2ip9t] * = - p 2 ^ + |^| 2 ^ . (2.19) 

Vfsi( r ) 



The operator V is 

P = -<9 2 + V (p) 



. « 2 2 4d(d + 2)r M + 2 - 4d(d + 3)r d+1 - (d - l) 2 
^O(P) = mV + -^ 16r ,- 1 ( rrf+1 _ 1) " L , (2-20) 



where r has to be expressed as a function of p using (2.16). 



The potential Vo(p) has the following behavior near the boundary and the tip 

m 2 , d(d+2) 

Vo(p) = 2-^ + 0(p 2 ), p^O, 

^O(P) = ~ 4( 1 )2 +Q ( X ) > ( 2 ' 21 ) 

The behavior at the boundary p = is of course the same as in pure AiSd+2. The behavior 
near the tip p = p* is in fact the correct behavior expected from a flat two dimensional 
space. Near the tip of the soliton the space becomes M 2 x with y = (p* — p) playing 
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the role of a radial variable and 6 playing the role of the polar angle. Indeed with the 
redefined field 

¥ y = -M= , (2.22) 



the operator V becomes, near y 
dimensions 



Vp* - p 

0, the zero angular momentum Laplacian in two 



y^O 



y 



•(V 2 



constant 



In fact the eigenvalues of the operator V which acts on \1/ 
V = -(V 2 ) 2 \ + V 1 (y) , 



V 1 (y)=V (y)+ A ., 



1 

4p 



(2.23) 



(2.24) 



are all positive. For d = 3 the proof is the following. Let us rewrite the potential V±(y) 
as follows 

V l (y) = (m 2 + ^)r 2 + V 2 (y) (2.25) 

where 



V 2 {y) 



1 + 3r 4 



+ 



(2.26) 



The term ^(y) is explicitly positive for all r. This may be seen as follows. The condition 
for positivity of V 2 (y) is 



•16 f 2 



- y(r) > 



(2.27) 



The inequality is saturated for y 
hand side becomes 



1 + 3r 4 

(r = 1). Furthermore the first derivative of the left 



3r 8 + 6r 4 _ x 
(1 + 3r 4 ) 3 / 2 



- 1 



VH^T L (1 + 3r 4 ) 3 / 2 1 (2 ' 28) 
This can be explicitly checked to be positive for all r > 1 (e.g. by squaring the expression). 
Therefore V^Q/) > for all r > 1. The first term in V±(y) in (2.25) is the asymptotic 
potential in AdS$ - when m 2 + > — | (which is the BF bound), this potential does not 
have any bound state. Since V^y) differs from this asymptotic potential by a positive 
function, the full potential V\(y) does not have any bound state. 

To look for a condensate in standard quantization, we need to find time independent 
solutions of the equation (2.19) which satisfy the boundary condition J = at p = and 
is regular at the tip p — p*. With these boundary conditions the operator V has a discrete 
and positive spectrum. This means that for sufficiently large p, the operator 



V EE V - p 2 , 

will have a negative eigenvalue. This is what we found numerically. 



(2.29) 
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At the critical value pi = /i c i the operator V has a zero eigenvalue, i.e. a zero mode 
which satisfies the appropriate boundary conditions both at the tip and at the boundary. 
This zero mode will play a key role in the following. 

Note that even though the operator V has negative eigenvalues in the condensed phase, 



the trivial solution does not become unstable. This is clear from (2.19) and from the fact 
the spectrum of V is positive, which shows that the frequencies of the solutions to the 
linearized equation are all real. 

Following the arguments of [16] it can be easily checked that the transition is standard 
mean field. This means that 

(O) J=0 ~ vVci - /"I , 
(0},=, cl ~ |J| 1/3 - (2-30) 

We expect that this transition extends to non-zero temperature, though we have not 
checked this explicitly. 

3 Quantum quench with a time dependent source 

We will now probe the quantum critical point by quantum quench with a time dependent 
homogeneous source J(t) for the dual operator O, with the chemical potential tuned to 
\i = /i c i. The function J(t) will be chosen to asymptote to constants at early and late 
times, e.g. 

J(t) = J tmh(vt) . (3.31) 

Note that we are using units with tq = 1. The system then crosses the equilibrium 
critical point at time t = 0. The idea is to start at some early time with initial conditions 
provided by the instantaneous solution and calculate the one point function (0(t)). In 
standard quantization this means that we impose a time dependent boundary condition as 



in (2.11) and calculate A(t). In alternative quantization the source should equal A(t). In 
this paper we discuss the problem in standard quantization : the treatment in alternative 
quantization is similar. 

3.1 Breakdown of adiabaticity 



With a J{t) of the form described above (e.g. (3.31)), one would expect that the initial 
time evolution is adiabatic for small v so long as Jo is not too small. As one approaches 
t — adiabaticity inevitably breaks down and the system gets excited. In this subsection 
we determine the manner in which this happens. 
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An adiabatic solution of (2.19) is of the form 



(p,t) = ¥°\p, J(t)) + e^ (1) (p,t) + eV 2 > + ■ ■ ■ , (3.32) 
where e is the adiabaticity parameter. The leading term is the instantaneous solution of 



(2.19), which is (using the definition (2.29)) 



P^ (0) + G(p)|^ (0) | 2 ^ (0) = , 



satisfying the required boundary condition. Here we have defined 

r 2— d 

G(p) 



(3.33) 



(3.34) 



From (2.30) we know that for a real J(t), this is real and has a form 
tt<°) ~ P a J(t) [1 + 0(p 2 )] + p x - a [J(t)] 1/3 [1 + 0{p 



where 



a 



d/2 . 



(3.35) 



(3.36) 



This follows from the equations (2.17), (2.18) and (2.11). The adiabatic expansion now 



proceeds by replacing d t — > ed t in (2.19) substituting (3.32) and equating terms order by 



order in e. The n-th order contribution ty( n > satisfies a linear, inhomogeneous ordinary 
differential equation with a source term which depends on the previous order solution 
i|r(n-i)_ lowest order we have the following equations for the real and imaginary parts 
of 



[V + 3G(p)(¥ 0) ) 2 ] (Re^ (1) ) = 0, 
[V + G(p)(¥ 0) ) 2 ] (lm¥ 1] ) = 2p,d t ¥ 0) 



(3.37) 



Note that in these equations the time dependence of J(t) should be ignored. The full 
function ^ must satisfy the boundary condition lim [p~ a ty(p,t)] = J(t). This means that 



the adiabatic corrections must start with the subleading terms, \E ,( - 1 ^ 



P 



l-a 



as p — > 



and has to be regular as p — > p*. These provide the boundary conditions for solving the 



equations (3.37). Consider first the equation for Im^ 1 ). Since the time dependence of 



\I/( ) is entirely through J(t) the solution may be written as 

Im*«(p,t) = 2^ J{t) J dp>g(p,p') ( p ',J( t )) 



(3.38) 
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where G(p,p') is the Green's function for the operator V + G{p)(^> 



(0)\2 



s(p, p') = wkm ^ ' p < ' ( 3 - 39 ) 

where and ip2 are solutions of the homogeneous equation [T> + G(p)(\l/^ ^) 2 ] ipi^ = 
which satisfy the appropriate boundary conditions at the tip p = p* and at the boundary 
p = respectively. The Wronskian W(ipi,i/)2) for this operator is clearly constant and is 
conveniently evaluated near the tip. Near p = p* these solutions behave as 

V>i ~ CVP* - P , ^2 ~ - P + BtJp~^p log(p* - p) , (3.41) 

where A, B, C are constants which depends on J(t) Thus the Wronskian is 

W{il>iM = -BC. (3.42) 

As noted in the previous section, the operator T> has a zero mode, i.e. [V + G(p)(^^) 2 ~\ 
has a zero mode when = 0, i.e. exactly at the equilibrium critical point. Thus, at 
this point we must have B = 0. This is why the first adiabatic correction lm^^(p, t) 
diverges at this point. For small J(t) we can use perturbation theory to estimate the 
value of B. For small J the zeroth order solution behaves as J 1//3 (the first term 



in (3.35) is subdominant). This is explicit to all orders in the expansion of the solution 
around the boundary. However, this is also justified by the results of the next section 
where we show that in the critical region the dynamics is dominated by a zero mode. 
The coefficient of the zero mode can be seen to be proportional to J 1//3 using a regularity 
argument similar to that in [16] so that the additional term in the operator behaves as 
G(p)(^(°)) 2 ~ [J(t)] 2/3 . This yields B ~ J 2/3 as well. Thus the Green's function which 



appears in (3.38) behaves as J 2 / 3 so that the correction behaves as 

(3.43) 



J 2 / 3 dJ(t) J 4 / 3 



The same argument shows that Re^ 1 ) = 0, so that \^^\ ~ ~jt/3 as well. Therefore 
adiabaticity breaks when 

|^(l) | _ |^(0) | ^ _ J5/3 _ ( 3>44 ) 



Note that in the equations (3.37) the time is simply a parameter 
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For sources which vanish linearly at t — 0, i.e. J(t) ~ ft (e.g. of the form (3.31)) this 
means that if the source is turned on at some early time, adiabaticity breaks at a time 



tadia ~ V~ 2/5 . (3.45) 

while at this time the value of the order parameter (O) is 

{0^)) ~ [J(tad»)] 1/3 = b^a] 1/3 ~ ^ • (3.46) 

With the usual adiabatic-diabatic assumption, these exponents lead to Kibble-Zurek scal- 
ing for a dynamical critical exponent z = 2, even though the underlying dynamics is 
relativistic and non-dissipative. From the above analysis it is clear that this happened 
because the leading adiabatic correction is provided by the chemical potential term, which 
multiplies a first order time derivative of the bulk field. 



3.2 Critical dynamics of the order parameter 

The breakdown of adiabaticity means that an expansion in time derivatives fail. In this 
subsection we show, following closely the treatment of 15 , that we now have a different 
small v expansion in fractional powers of v during the period when the sources passes 
through zero. This will lead to a scaling form of the order parameter in the critical 
region. In the following we will demonstrate this for the case where J(t) ~ vt near t ~ 0. 
However the treatment can be easily generalized to a J(t) ~ (v t) n for any integer n. 
To establish this, it is convenient to separate out the source term in the field ^(p,t), 



y(p,t) = p a J(t) + y s (p,t) 



a 



A_ - al/2 



(3.47) 



where we have used the relation (2.18) and the fact that near the boundary p ~ 1/r. The 
equation of motion (2.19) then becomes 

- d^ s + 2ipd t ^ s = (Vp a )J(t) + Vm s + G{p) [p 3a [J(t)f + p 2a [J(t)] 2 (2^ s + % 

+G(p) [p a J(t)(2\y s \ 2 + y 2 s ) + \y s \ 2 V s ] 



+p a [dfJ-2ipd t J] . 



(3.48) 



This separation is useful because we know that in the presence of a constant source 
J(t) = J, the static solution has the asymptotic form 



p l-a pi/3 + 0(/j2) ] + Jf) a + 2 ^ + Q ^ 



(3.49) 



which follows from (2.30). 
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The scaling relations (3.45) and (3.46) suggest that we perform the following rescaling 
of the time and the field 



t = v'^rj 



= v 1/5 X ■ 



(3.50) 



In the critical region we can now use J(t) = vt = v 3 ^ 5 i] and rewrite (3.48) as an expansion 
in powers of v 2 ^ 5 , 



V X = v 2 ^ [2ipd vX - G{p)\x\\ ~ v(Vp a )} + 0(i 



,4/5^ 



(3.51) 



As noted above, because of the boundary condition at p = and the regularity condition 
at p = p* the spectrum of V is discrete. Let (p n be the orthonormal set of eigenfunctions 
of the operator T> 

Vip n (p) = X n (fn(p) , n = 0, 1, • • • , (3.52) 

with A = 0. (po(p) is the zero mode which we discussed earlier. Since p has been tuned 
to be equal to p cl , all the higher eigenvalues are positive. 
We now expand 

X(P,V) = ^2xn(r))<Pn(p) , (3.53) 



and rewrite the equation (3.51) in terms of the modes Xn(v) 



2/5 



C ni n 2 n 3 Xn 3 Xn2Xni 4" ^TnV 



rtira 2 "-3 



+ 0(v 



4/5 \ 



(3.54) 



where we have defined 



yJn 



f" 



dp^{ P ){Vp a ) , 

0)<Pn 2 (p (p)G(p) . 



(3.55) 



It is clear from (|3.54|) that the zero mode part of the bulk field dominates the dynamics 

(3.56) 



in the critical region. In fact for small v a solution is of the form 

Xn(v)=5noUv)+V 2/ % + 0(v 4 / 5 ) . 

The zero mode satisfies az = 2 Landau-Ginsburg equation 

- 2ipd^ + e o o |£or£o + JoV = 0. 
Reverting back to the original variables we therefore have 



(3.57) 



(3.58) 
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which implies a Kibble-Zurek scaling for the order parameter with z = 2 



(0(t,v)) =v 1/5 (0{v 2/5 t,l)) . (3.59) 



Note that the effective Landau-Ginsburg equation (3.57) is not dissipative because the 
first order time derivative is multiplied by a purely imaginary constant. In fact, in the 
absence of a source term the quantity ^(|£o| 2 ) 2 is independent of time. 

Beyond the critical region, we cannot use the approximation J(t) ~ vt and there is 
no useful simplification in terms of the zero mode. However, the boundary conditions at 
the tip are perfectly reflecting boundary conditions (as appropriate for the origin of polar 
coordinates in two dimensions) so that there is a conserved energy in the problem. This 
is in contrast to a black hole background where there is a net ingoing flux at the horizon 
causing the system to be dissipative. Indeed in the quench problem considered in [15] 
arguments similar to those used in this section also led to an effective Landau-Ginsburg 
dynamics with z = 2, but which is dissipative. 

In the appendix we analyze a Landau-Ginsburg toy model motivated by the results of 
this section. 



4 Numerical results 

In this section we summarize our numerical results. We have solved the bulk equation 
of motion numerically for d = 3. The results for different values of m 2 are similar. We 
present detailed results for m 2 = — 15/4. In this case the critical value of the chemical 
potential is \i c \q ~ 1.88. 



We discretize the partial differential equations (PDEs) (2.19) (written in the y-coordinate) 
in a radial Chebyshev grid to study the numerical problem. Once discretized in radial 
direction, the PDEs become a series of ordinary differential equations (ODEs) in the 
temporal variable. The resulting ODEs are solved with a standard ODE solver (e.g. 



CVODE). The time dependence is chosen to be of the form dependent source as in (3.31 ). 
In principle one may study with any kind of time dependent source. 

We will consider the problem in two regimes. The first is the "slow" regime where we 
expect our analytic arguments to be accurate, the other is a "fast" regime where there 
is no adiabatic region whatsoever. In the slow regime we will try to zoom on the scaling 
region around the phase transition. In the fast regime we will find a large deviation from 
the adiabatic behavior and possible chaotic behavior. 
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4.1 Slow regime 

Since our main interest is quench through the critical point, we concentrate mainly near 
the phase transition. We choose \iq = /i c ig(~ 1.88), so that the system is critical in the 



absence of any source. In the presence of a time dependent source of the form (3.31) 
we calculate the bulk field ^(t) and extract from this the value of (0(t)) of the dual 
field theory. A typical plot of the real part of (0(t)) for slow quench through the phase 
transition is presented in Figure |2j 



Re<0(t)> 





Figure 2: The plot of Re (0(t)) with v = 0.02. 

Clearly the late time behavior is oscillatory, reflecting the fact that we are dealing 
with a closed and non-dissipative system. 

We then zoom on the critical region near t = for various value of v to look for any 
scaling behavior. One way to look for this is to consider the behavior of (O) at t — 0. 
Equation ( |3.59[ ) then predicts a scaling behavior (O(0)) ~ v 1 ^ 5 . 

Figure [3] shows a plot of log(Re (O(0))) for different v. We fit the data points with 
a function f(x) = A + Bx + C/x, where x is log(f). Here we kept a sublinear (0(l/x)) 
term to understand how the fit function approaches a linear regime. From our analytic 
argument we expect B = 1/5. A fit of the numerical results yields f(x) = 0.794 — 
0.490/x + 0.206x. Changing the number of fit points and range changes the values of fit 
parameters a bit, however we always get a value of B which is close to 1/5 with only a 
few percentage deviation. The imaginary part (Im (O(0))) also satisfies the same scaling. 

4.2 Fast regime 

In the fast regime we see a large deviation from the adiabatic behavior, as shown in Figure 
H 
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Log(Re<O(0)>) 




Figure 3: The plot of log(Re (O(0))) vs log(v). We also plotted the closest fit (see text). 

Re<0(t)> 




Figure 4: Plot of Re (0(t)) showing chaotic behavior with a large value of v = 2.0 and 



The motion in this regime becomes possibly chaotic. Here we have a system with a 
conserved energy. Once we put some energy in the system, the non-linearity possibly takes 
the system over the whole phase space (Arnold diffusion). It is expected that if we wait 



long enough the probe approximation actually breaks down 26 and we have to consider 



the fully backreacted problem. We plan to attack this problem in the near future. 
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A Adiabatic and scaling analysis of a toy model 

In this appendix, we consider a (0 + 1) -dimensional toy model which follows the equation 

2i/i0+ (m 2 -/i 2 )0 + 0|0| 2 = J(t) . (A.60) 

The function J(t) asymptotes to constants at early and late times and passes through 
zero in a linear fashion at some intermediate time, e.g. 

J(t) = J tanh«) . (A.61) 

A.l Adiabaticity 

We first derive conditions for breakdown of adiabaticity near the critical point m 2 = /i 2 
and J(t) = 0. We carry out adiabatic expansion as following: 

d t ^ed t , 0-+0 o (J(t)) + e 1 (t) + ... , (A.62) 

where <f>o(J(t)) is the (real) adiabatic solution given by 

MJ(t)) = [J(t)] 1/3 ■ (A.63) 

The solution to 0(e 2 ) is , 

• 2 

^ = MAt)] + ei |0 O + e 2 -^ (s^A - ^) - ^0 O 2 ) + 0( e 3 ) . (A.64) 
where The breakdown of adiabaticity happens when, 



(A.65) 



For J(t) = tanh(t>t) ~ vt the above two equations translate into, 

H ~ . (A.67) 
Thus if \i is of 0(1) then the above equations give us, 

t ~ v 2/b . (A.68) 
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A. 2 Scaling behavior 

Now sitting at the critical point we study the behavior of the scaling solution with fi = 
0(1). From the adiabatic analysis we expect scaled time, i = v 2 ^t. 

We write the field as x + an d the source as Jr + iJ\ m , where both Jr and J\ m go as 
vt. To find the scaling exponents we extract the v dependencies as, 

t = v a t, X = v? X , Z = v^. (A.69) 
Consistency of the equations demand, 

« = -?. P = l, 1=\- (A.70) 
5 5 5 

This determines the scaling behavior of the field (p at m 2 = fi 2 and with \i of 0(1) as, 

cj)(t, v ) = v 1/5 (j)(v 2/5 t,l) . (A.71) 

This agrees with our expectation from adiabatic analysis, and also has been confirmed 
numerically in Section 4. 



A. 3 Late time behavior 

At late times the source J(t) can be treated to be a constant and the solution can be 
obtained by perturbing around the static solution, 4> sta tic = J 1 ^ 3 - It is then straightforward 

t2/3 

to see that the solution 4>{t) is oscillatory with a frequency uj — , 



References 

[1] For reviews and references see S. Mondal, D. Sen and K. Sengupta, arXiv:0908.2922; 
J. Dziarmaga, arXiv:0912.4034; A. Polkovnikov, K. Sengupta,A. Silva and M. Ven- 
galattore, dynamics of closed interacting quantum systems" , |arXiv: 1007.5331 



[2] P. C alabrese and J. L. Cardy, J. Stat. Mech. 0504 (2005) P010 |arXiv:cond-| 
mat/0503393 1; P. Calabrese and J. L. Cardy, Phys. Rev. Lett. 96 (2006) 136801 
|arXi v : cond-mat / 0601225] 



[3] P. Calabrese and J. Cardy, |arXiv:0704.1880 [cond-mat. stat-mech]]; S. Sotiriadis and 
J. Cardy, J. Stat. Mech. (2008) P11003T |arXiv:0808.011~6~| [cond -mat.stat-mechll; 
S. Sotiriadis, P. Calabrese and J. Cardy, EPL 8 7 (20097^ 0002, | arXiv:0903. 0895 
[cond-mat. stat-mech]]; S. Sotiriadis and J. Cardy, [arXiv:1002.0167 [quant-ph] . 



[4] J. M. Malda cena, Adv. Theor. Mat h. Phys. 2 (1998) 231 [Int. J. Theor. Phys. 38 
(1999) 1113] larXiv:hep-th/9711200] ; 



20 



[5] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998) 
|arXiv:hep-th/9802109] . 

[6] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998) |arXiv:hep-th/9802150| . 

[7] O. Aharony, S. S. Gubs er, J. M. M aldacena, H. Ooguri and Y. Oz, Phys. Rept. 323, 
183 (2000) |arXiv:hep-th/9905111| . 

[8] R. A. Janik and R. B. Peschanski, Phys. Rev. D 74, 046 007 (2006) |arXiv:hep-| 
|th/0606149| ; R. A. Janik, Phys. Rev. Lett. 98, 022302 (2007) |arXiv:hep-t h/0610144]7 
P. M. Chesfer and L. G. Yaffe, supersymmetric Phys. Rev. Lett. 102, 211601 (2009) 
|arXiv:0812. 2053 [hep-th]]; P. M. Chesler, L. G. Yaffe, dynami cs in N = 4 supersym- 
metric Yang-Mills theory," Phys. Rev. D82, 026006 (2010). |arXiv:0906.442~6] [hep- 
th]]; S. Bhattacharyya and S. Minwalla, JHEP 0909 (2009) 034 |arXiv: 0904.0464] 
[hep-th]]. D. G arfinkle, L. A. Pando Zayas, Phys. Rev. D84, 066006 (2011). 
larXiv:1106.2 339 [hep-th]]. 

[9] J. Abajo-Arrastia, J. Aparicio, E. L opez, JHEP 1 011, 149 (2010). |arXiv:1006.4090| 
[hep-th]] ;H. Ebrahim, M. Headrick, |arXiv:1010.5443 [hep-th]]; V. Balasubramanian, 
A. Bernamonti, J. de Boer, N. Copland, B. Craps, E. Keski-Vakkuri, B. Muller, 
A. Schafer et aL, Phys. Rev. Lett. 106, 191601 (2011). |arXiv:1012.4753| [hep- 
th]]; V. Balasubramanian, A. Bernamonti, J. de Boer, N. Copland, B. Craps, 
Keski-Vakkuri, B. Muller, A. Schafer et al, Phys. Rev. D84, 026010 (2011). 



arXiv:1103.2683 [hep-th]]; E. Caceres and A. Kundu, JHEP 1209, 055 (2012) 
arXiv:12 05.2354 [hep-th]]; D. Galante and M. Schvellinger, JHEP 1207, 096 (2012) 
arXiv:1205.1548 [hep-th]]. 



[10] T. Alba sh, C. V. Johnson, Electromagnetic Quenches," New J. Phys. 13, 045017 
(2011). |arXiv:1008.3027 [hep-th]]; T. Takayanagi, T. Ugajin, Coarse-Graining," 
JHE P 1011, 054 (2010). [arXiv:1008.3439] [hep-th]]; C. T. Asplund, S. G. Avery, 
larXiv:1108.25l0| [hep-th]]. 

[11] S. R. Das, T. Nishioka, T. Takayanagi, JHEP 1007, 071 (2010). |arX iv:1005.3348| 
[hep-th]]. 

[12] A. O'Bannon, |arXiv:0808.1115| [hep-th]. K. Hashimoto, N. Iizuka, T. Oka, 
|arXiv:1012.4463| |hep-th ||; K. -Y Kim, J. P. Shock, J. Tarrio, JHEP 1106, 017 
(2011). |arXiv: 1103.45811 [hep-th]]; S. Prem Kumar, Phys Rev. D 84, 026003 (2011). 
jarXiv: 1104. 1405 [hep-th]]. S. Janisz ewski, A. Karch [arXiv:1106.40lo"l [hep-th]]; 
C. Hoyos, T. Nishioka, A. O'Bannon, |arXiv: 1106.4030] [hep-th]]. 

[13] T. W. B. Kibble, J. Phys. A A 9, 1387 (1976); W. H. Zurek, Nature 317, 505 (1985). 

[14] A. Chandran, A. Erez, S. Gubser and S. Sondhi, Phys. Rev. B 86, 064304 (2012). 

[15] P. Basu and S. R. Das, JHEP 1201, 103 (2012) [arXiv: 1109.3909] [hep-th]]. 

[16] N. Iqbal, H. Liu, M. Mezei, Q. Si, Phys. Rev. D82, 045002 (2010). | arXiv:1003.001~0| 
[hep-th]]. 

[17] S. S. Gubser, Phys. Rev. D78, 065034 (2008). |arXiv: 080 1.2977] [hep-th]]. 



21 



[18] S. A. Hartnoll, C. P. Herzog, G. T. Horowitz, JHEP 0812, 015 (2008). 
|arXiv:0810.l "563|[hep-th11; S. A. Hartnoll , C. P. Herzog, G. T. Horowitz, Phys. Rev. 
Lett. 101, 031601 (2008). [arXiv: 0803. 3295] [hep-th]]. 

[19] P. Basu, A. M ukherjee, H. -H. Shieh, Phys. Rev. D79, 045010 (2009). 
| arXiv:0809.4494| [hep-th]]; 

C. P. Herzog, P . K. Kovtun, D. T. Son, Phys. Rev. D79, 066002 (2009). 
larXiv:0809.4870 [hep-th]]. 



[20] S. A. Hartnoll, Class. Quant. Grav. 26 (2009) 224002 |arXiv:0903.3246 [hep-th]]. . 



[21] D. Arean, P. Basu, C. Krishnan, JHEP 1010, 006 (2010). |arXiv:1006.5T65] [hep-th]]; 

[22] C. P. Herzog, J. Phys. A 42, 343001 (2009) [arXiv:0904.1975| [hep-th]]; 

[23] For review and further references see T. Faulkner, N. Iqbal, H. Liu, J. McGreevy, 
D. Vegh, |arXiv: 1101.0597 [hep-th]]; 



[24] T. Nishioka, S. Ryu and T. Takayanagi, JHEP 1003, 131 (2010) [arXiv: 091 1 .0962 
[hep-th]]. 

[25] K. Murata, S. Kinoshita and N. Tanahashi, JHEP 1007, 050 (2010) |arXiv:1005.0~6~33| 
[hep -th] ]; M. J. Bhaseen, J. P. Gauntlett, B. D. Simons, J. Sonner and T. Wiseman, 
|arXiv: 1207.4194] [hep-th] ; 



[26] P. Bizon and A. Rostworowski, Phys. Rev. Lett. 107, 031102 (2011) |arXiv:1104.3702 



[gr-qc]]; O. J. C. Dias, G . T. H orowitz, D. Marolf and J. E. Santos, Class. Quant. 
Grav. 29, 235019 (2012) |arXiv: 1208.5772] [gr-qc]]. 



22 



